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Abstract 

Recent results from the HESS gamma ray telescope have shown the presence of 
both a diffuse, extended, flux of gamma rays above ~0.4 TeV and discrete sources 
in and near the Galactic Centre. Here, we put forward a possible explanation in 
terms of the diffusion of cosmic ray protons from a succession of supernova remnants 
( SNR ) in the SgrA* region of the Galaxy plus a contribution from SNR in the 
rest of the Galactic Centre Region, to be called the Galactic Centre Ridge ( GCR ). 
Protons are favoured over electrons because the mG magnetic fields in the Region 
will attenuate energetic electrons severely. 

Prominent features are the need for 'anomalous diffusion' of the protons in the 
whole region and the adoption of low efficiency for SNR acceleration in the high 
density regions. The latter is related by us to the well-known low 'cosmic ray 
gradient' in the Galaxy. 

A corroborating feature is the close correlation of inferred cosmic ray intensity 
with the smoothed intensity of 5 GHz radio radiation. We attribute this to the 
presence of the SNR in the GCR. 

1 The Galactic Centre 

It is well known that the Galactic Centre Region ( GCR ) has remarkable properties, and 
it is not surprising that the Cosmic Ray ( CR ) view has singular interest. 

Concerning molecular gas, some 10% of such gas in the Galaxy is to be found in the 
central 200 pc radius ( Morris and Serabyn, 1996, to be denoted I ). The magnitude is 
quoted as (5 - 10) ■ lO'^ Mq (I), with probably 4.4 ■ lO^M© in the region of concern here, 
r < 200 pc (Tsuboi et al., 1999). The mean density is very high, (n) > lO'' cm~^, in 
many of the clouds and the filling factor similarly high: ~ 0.1. The mean temperature is 
elevated, being typically 70K ( in the range 30K - 200K ) (I). Surrounding the molecular 
clouds is a very hot plasma of 'temperature' 10-15 keV (I), a value which is so high 
as to need the presence in the past of very energetic activity, such as a giant Galactic 
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Centre explosion and/or the explosion of many SN within the central parsec over the past 
10^- 10^ years (I). 

Many measurements have indicated the existence of a magnetic field of several mG, 
with both tangled and smooth components. A prominent feature is the presence of field 
'tubes' nearly perpendicular to the Galactic Plane (I). 

SgrA* comprises a compact central object containing a radio source with a fiat spec- 
trum. It is almost certain that a black hole of mass (2 — 3) x 10^ Mq is involved (I). 

Of many other unusual features associated with the GCR, mention should be made 
of the high velocity winds ( 500 - 1000 kms~^ ) ( Breitschwerdt et al.,2002 ; Volk and 
Zirakashvih,2004 ) associated particularly with the region SgrB2 ( eg Sunyaev et al., 
1993 ). 

Particular remarks should be directed to the 'energetics' of the GGR; we are mindful 
of the fact that, locally, the energy densities of GR, starhght, magnetic fields and gas 
motion are all equal at ~ Q.beVcmr^ ( Wolfendale, 1983 ). In the GGR we estimate, 
from the published data given in (I), that the energy densities of plasma, the far- infra-red 
radiation, magnetic fields and gas motion are all several thousand eVcm~^ ( indeed, the 
magnetic energy density will be higher still over limited volumes ). 

The features listed above all lead to particular interest in the gamma ray signal from 
the general region of the Galactic Gentre. Early work by SAS II ( Thompson et al.,1976 ), 
GOSB ( Mayer-Hasselwandler et al.,1982 ) and GGRO ( Mayer-Hasselwandler et al.,1998; 
Hartmann et al.,1999 ) in the 100 MeV region showed that the GR intensity there was 
probably not very different from that locally, but the angular resolution was too poor 
for detailed study. More recently, GANGAROO ( Tsuchiya et al.,2004 ), VERITAS 
( Kosack et al.,2004 ), HESS ( Aharonian et al.,2004 ) and MAGIG ( Albert et al.,2006 ) 
have observed the GC in sub-TeV and TeV gamma rays. Now, HESS, with its superior 
angular resolution ( ~ 0.1° ) coupled with its higher energy range ( from hundreds of 
GeV to tens of TeV ), offers the possibility of another window on this dramatic region 
( Aharonian et al., 2006 ) and we shall base our analysis mainly on this latest work. 
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2 The cosmic ray aspect 



Over the years we have examined the imphcations of cosmic ray acceleration in supernova 
remnants (SNR) for a wide variety of CR phenomena. A strong case has been made for 
such acceleration providing the bulk of CR up to 10^^ eV (eg Erlykin and Wolfendale, 
2001) and perhaps - for a special class of SNR - beyond (Biermann, 1993; Sveshnikova, 
2003; Erlykin and Wolfendale, 2005 ). 

Gamma ray astronomy provides a proxy indicator of CR in distant parts of the Galaxy 
by way of gamma ray production in CR - interstellar gas collisions and CR interaction with 
magnetic and photon fields. A problem is that the nature of the CR (nuclei or electrons) 
is not known, a priori, and also there are often difficulties with the determination of 
the necessary column density of gas. Nevertheless, in the absence of better ways of 
studying the origin of CR, we proceed, noting particularly that high energy electrons - 
those necessary for TeV gamma rays - will be severely inhibited by the niG magnetic fields, 
whether they be primary or secondary ( in any event, the fiux of secondary electrons well 
away from the source will be very small ). 

HESS has given contours of the TeV gamma ray emission from the exciting Galactic 
Centre Region ( GCR ), specifically, that bounded by -2° < I < +2°, -1.2° < h < 1° 
( http: / / www.mpi-hd.mpg.de/hfm/HESS ) . The 'map' is characterized by 'point source' 
emission from the 'compound SNR' 6*0.9 + 0.1 and from SgrA*, very close to the nominal 
GC. 

3 The basic data 
3.1 The gamma ray map 

Following the analysis of the HESS group, we consider the region: |/| < 1.6° and |&| < 0.3°. 
Figure 1 shows the profile of the 'intensity' (actually the counts) as a function of longitude 
( Aharonian et al.,2006 ). The result of subtracting the two point sources {G 0.9 + 0.1 
and Sgr A*) is indicated. The workers used the known point spread function (PSF) for 
this subtraction. The background has been also subtracted. 

The authors point out that the diffuse gamma-ray emission diminishes with longitude 
in the studied region and fades away at |/| > 1.2° inspite of the fact that there is a 
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substantial amount of gas at least at positive longitudes of I > 1.2°. They explain this 
profile assuming that gamma rays are produced in CR - gas collisions and CR have a 
gaussian distribution around the GC with a best fit width of 0.8°. 

In what follows we draw the physical scenario of the observed features. We examine 
the possibility that our SNR model has validity here, too, viz that SNR are the sources 
and that the CR propagate in a diffusive manner. Specifically, we assume that many SNR 
occurred near the GC during the last 10^ — 10^ years, ( and extending back, perhaps, even 
longer ) as already suggested. As in the HESS group scenario the gamma rays came from 
the interactions of CR protons from SNR with the ambient high density gas, essentially 
all the gamma rays from SgrA* coming from the burst of SNR, the remainder ( from the 
GCR ) coming either from the burst of SNR or the 'conventional SN' in the GCR ( which 
would surely be expected in view of the considerable amount of stellar activity ). 




1 -1 

Galactic Longitude (degrees) 

Fi gure 1: The HESS results on gamma rays from the Galactic Centre Region ( Aharonian et al., 2006 ). 
The histogram gives the profile of gamma ray intensity above 0.38 TeV for |6| < 0.2°. Two 'point' sources 
have been removed, as indicated by dashed lines. The solid line is the column density of molecular 
hydrogen between the same latitude limits normalised to the same total area. The dashed line shows the 
gamma-ray flux expected if the CR density distribution can be described by a Gaussian centred at Z = 0° 
and with rms 0.8° ( following the analysis of the HESS group ). 
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3.2 The distribution of target gas and the inferred CR intensity 

Many measurements have been made of the distribution of gas in the GC region, most 
notably using CO (eg Bania et al., 1977, Oka et al., 1998) and CS ( Tsuboi et al., 1999). 
An apparent feature, common to all analyses, is the very high density of gas (^^2) in the 
region. As the last authors have remarked, there are peculiar structures here which are 
'presumably related to the unique activity in the GC region'. The HESS group used the 
CS data to give the consequent column density of molecular hydrogen and we have done 
the same, although we have applied a correction for the loss of lower density molecular 
gas from the work of Dame et al.,2001 for the region \l\ < —1°. 

At the beginning we tried to find the general characteristics of CR in the GC region. 
Unlike the HESS group's assumptions about the gaussian distribution of CR around the 
GC we adopted a uniform distribution in this area and started with the local value of 
the CR intensity. Using our emissivity of the yield of gamma rays per hydrogen atom 
for the local CR spectrum (Erlykin and Wolfendalc, 2003) we have derived R(l) - the 
ratio of the observed/expected gamma ray intensity with the result shown in Figure 2. If 
the observed gamma rays originate indeed from CR - gas collisions the gas distribution 
is eliminated from R{1) and this ratio gives the actual CR longitudinal profile compared 
with the assumed uniform distribution. The possible contribution of electrons is ignored 
in this assumption following the arguments in §2. 

The mean value of {R{1)) averaged over the interval of -1.4° < / < 1.4° is 2.44 ± 0.41 
for the total flux including the central source. Without the central source it is {R{1)) = 
1.83 ± 0.15. It is with the interpretation of this plot that we shall be mainly concerned. 

3.3 The radio map 

Of likely relevance is the radio map for the same region. Figure 2 shows the average 
intensity over \b\ < 0.3" derived by us from the 5 GHz data of Altenhof et al., 1979, that, 
for 10 GHz is very similar ( Handa et al.,1987 ). Its relevance will be considered briefly 
here and in more detail in §5.3. 

Inspection of Figure 2 shows that there is a correlation between the excess CR intensity 
and the radio emission intensity and this is quantifled in Figure 3. The slope of the straight 
line flt is 0.96±0.11 and the corrrelation coefficient is 0.85; there is thus strong evidence 
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full line: ratio for total gamma rays 
dashed line: ratio for diffuse gamma rays 
dotted line: radio ( 5 GHz ), arb. units 



1 .5 




-1 .5 



Figure 2: Ratio of observed to expected gamma ray intensity for E-^ > 0.38 TeV, the expectation being 
that the CR intensity is the same as that locally; the ratio is thus that of the proton intensity to that 
locally: R{1). The full line shows the ratio for the total gamma-ray intensity, the dashed line - the same, 
but with the central source subtracted, the dotted line - the profile of the radio intensity in arbitrary 
units. The contribution from CR electrons is ignored following the arguments in §2. Inspection of the 
basic data in Figure 1 shows that most of the peaks are significant at the 2-3 standard deviation level. 



for linearity. 

Although some of the radio photons come from discrete sources ( 'young' pulsars, etc ) 
and SNR at the 'centers' of the condensed contours many come from surrounding regions 
and are due to CR electrons ( from old pulsars, SNR etc ) undergoing synchrotron radi- 
ation in the strong magnetic fields in the region, together with a 'thermal contribution', 
which in fact is greater than the non-thermal ( ie from synchrotron radiation ) at this 
frequency. 

In fact, the thermal contribution, too, is ( statistically ) correlated with the sites of 
SNR from reasons of energetics. 

Of importance is the fact that the positions of the peak intensities in the maps at 5 



GHz, 10 GHz, 60fi ( IRAS ), 10 GHZ, 5 GHZ and 2.7 GHz ( Effelsberg: |http: / /www. mpifr.de/old_mpifr/ 
are all coincident with the 'peaks' in the CR plot ( Figure 2 ). They are not coincident 
with the peaks in the column density of gas ( Figure 1 ). 

The likelihood of the radio flux being an indicator of the distribution of SNR ( 'present' 
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Figure 3: Correlation of the CR intensity ( as distinct from the gamma-ray intensity ) with the integrated 
radio intensity ( at 5 GHz ) for |&| < 3°. The results relate to the GCR. 

and past ) in the GCR is enhanced by the hnear size distribution of the excesses in Figure 
2, as will be demonstrated later. 

4 Analysis of the data 
4.1 General Remarks 

We start with Figure 1. It is interesting, and perhaps very significant, to note the large 
discrepancy in the region between I — +1° and I — +2° between the column density of 
gas and the measured gamma ray intensity - a discrepancy that would not exist if the CR 
intensity were constant. In other words, there is a significant fraction of the gas in the 
GCR ( some 18% ) that is under-populated by CR. It is this mass of gas, beyond I — 1°, 
that makes the gas non-symmetric about the Galactic Centre. The question of symmetry 
in the various GCR properties is taken up again later. 
Moving to Figure 2, a number of remarks can be made. 

(i) The likelihood of systematic errors in the inferred column density of molecular 
hydrogen means that the absolute values of the enhancement in CR intensity ratio 
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( denoted R{1) ) are uncertain. However, the shape of the longitude-dependence 
should be reasonable. 

(ii) There should be a trend of diminishing R (/) with increasing |/| if our contention 
about diffusion of CR from sources in SgrA* is correct. 

(iii) Some measure of correlation of R (/) with the integrated radio intensity should occur 
if, as has been remarked, there arc SNR in the GCR itself. An alternative way of 
approaching this problem is in terms of CR diffusion - in high B regions there will 
be a smaller diffusion coefficient because of increased turbulence. 

4.2 Overall CR excess 

Inspection of Figure 1 shows that the gamma ray excess falls down slightly when receding 
from the Galactic Centre, ie with increasing |/| and, correspondingly, the CR intensity also 
falls. Certainly the GCR is a singular region with unique properties. Any extrapolation 
of its characteristics to the wider longitude range should be taken with care. However, 
taking into account the large errors in the observed intensity of the diffuse gamma rays 
and uncertainty in the column density of the gas we cannot rule out a weaker longitude 
dependence of the CR intensity in the wider longitude range of the Inner Galaxy, than is 
indicated in Figure 2. 

The analysis of the diffuse gamma-ray profile at lower GeV energies obtained with 
the SASH satellite indicated a weak radial CR gradient in the Inner Galaxy, ( Issa and 
Wolfendale,1981 ), as did the CGRO satellite ( Strong and Mattox,1996 ). The CR 
intensity in the Inner Galaxy does not exceed the local value by a factor of more than 
1.2. If our latter value of the excess equal to 1.83±0.15 can be extrapolated into a wider 
longitude region than |/| < 1.6°, it might mean that the CR gradient for sub-TeV and 
TeV energies is stronger than for GeV energies, as has been predicted by us ( Erlykin and 
Wolfendale, 2002 ). 

4.3 The GC Region 

Our model, to be tested, comprises SN which have exploded in the central region ( to- 
gether, probably, with the 'normal quota' of SNR associated with the gas in general ), 
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providing protons which diffuse away and permeate the molecular gas, in which they pro- 
duce secondary pions and thereby gamma rays. It is supposed that recent SN provide the 
peak at the position of Sgr A* itself (to be discussed later) . 

In Figure 4 we show the gamma-ray profile expected from a single SN which exploded 
either 10^ or 10'^ years ago. Calculations have been made for both 'normal' and 'anoma- 
lous' diffusion ( the distinction relates to the nature of the 'scattering centres' in the 
interstellar medium ). Briefly, the difference between them can be reduced to two basic 
features: 

(i) the diffusion radius Rd{t) for anomalous diffusion depends linearly on time as Rd{t) oc t 
while for normal diffusion Rd{t) oc y/i. Both radii are normalized to 1 kpc at the time 
equal to the mean life time of CR particles for the particular energy; 

(ii) the lateral distribution function of the CR density p(r) for anomalous diffusion is not 
gaussian as for normal diffusion, but has a more complicated shape: p(r) oc (-^)^)~^. 
The difference between these two diffusion modes is discussed in more detail in Lagutin,2001, 
Erlykin and Wolfendale, 2002, Erlykin, Lagutin and Wolfendale, 2003. 

The effective averaged gas density in our calculations for the region in question |/| < 
1.6°, \b\ < 0.3° is lOOcm"^ ( this being much less than the high values of ~ lO'^cmr'^ quoted 
because of the presence in our studied volume of a large volume of low density material 
and the 'filling factor', ie the fraction of space occupied by the high density molecular 
gas ). 

It is seen that there is no sharp peak in the centre even for the youngest SN with an 
age of 10^ years. The fall of the intensity with longitude in the wings is too strong to be 
compatible with the experiment. If we reduce the age below 10^ years the sharpness of 
the peak increases but the fall of the intensity in the wings increases too. Therefore, to 

get agreement with experiment both in the central peak and in the wings ( Figures 1 and 
2 ) we cannot use just a single SN and need the succession of SN explosions distributed 
in time. 

In Figure 5 we show the results of the calculation for a succession of SN in Sgr A*. 
Two sets of SN rates are taken: 10^ SN in 10^ years and 10^ SN in 10^ years with the 
uniform temporal distribution of SN explosions within these time intervals. Since the 
average rate in both cases is 1 SN in 100 years the diffusion approximation for CR from 
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normal diffusion: 



anomalous diffusion 



full line: age = 10 yeors 
dotted line: age = 10^ years 
dashed line: age — 10*^ years 



dash — dotted line: age ^ 10 years 
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Figure 4: Longitude dependence of the gamma ray intensity from a single source at longitude ^, in a 
medium of density lOOcm^^ at 6 = 0°. Two ages are considered, 10^ and 10^ years and two modes of 
diffusion: normal and anomalous. 



some young SN and large distances from the GC cannot be valid because the CR velocity 
cannot exceed the speed of light. We have introduced limitations for this effect and have 
found that the gamma-ray intensity a,t i — 2° decreased by ~ 1.6%. In all subsequent 
calculations we applied these limitaions. 

In Figure 5 the ( important ) sharp spike, which results mainly from SNR younger 
than about 10^ years, from which the CR do not diffuse very far, is a consequence of 
anomalous diffusion. Normal diffusion gives a much weaker spike ( see also Figure 4 ). 
In practice, in view of the sources being in SgrA* itself, where the density is very much 
higher than lOOcm"^, the peak will be even higher. The other, later SN will have given 
particles which have diffused out to permeate more of the molecular material in the GCR. 

The extent to which the idealised 'lateral distributions' in Figure 5 would be modified 
using the actual column density of gas ( rather than gas of constant density ), averaging 
over the galactic latitude interval \b\ < 0.3° and applying the correction for the finite 
angular resolution of the telescope of 0.1°, can be seen in Figure 6. The general trend is 
seen to be preserved. For the comparison we also provide the results of HESS observa- 
tions. The difference between the absolute values of calculated and observed intensities 
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Figure 5: The angular profile of the gamma ray intensity in the GCR predicted for our model with 
n = lOOcm"^ and either 10'^ SN in 10^ years at the position of SgrA*, placed at the longitude £ = 0°, or 
10^ SN in 10* years. Results are given for boh anomalous and normal diffusion. 

is discussed below. 

5 Interpretation of the results 

5.1 General Remarks 

Guidance as to the frequency of SN comes from remarks in a variety of works ( eg I ) that 
of order 10^ SN in 10^ years to 10^ SN over the past 10^ years would have been sufficient 
to provide the energy necessary for the very strong wind and other features visible in this 
unique region of the Galaxy. Production over a longer period seems unlikely although it 
must be said that there is some evidence for a bout of star formation between 3-10^ and 
7-10^ years age ( Krabbe et al., 1995 ). 

It must be admitted that the manner in which the particles diffuse in this region is 
debatable both by way of the diffusion coefficient to adopt and the manner of diffusion, 
viz 'normal' or 'anomalous'. We consider that, in view of the disturbed conditions in the 
region and highly non-uniform distribution of gas, the mode of diffusion in the very central 
region at least will be 'anomalous' (Erlykin, Lagutin and Wolfendale, 2003). Goncerning 
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COMPARISON OF CALCUU\TIONS WITH OBSERVATIONS 
(10^ SNR in 10® years, normal and anomalous diffusion, |b| < 0.3° ) 

— I 1 1 1 1 1 

full line: HESS observation 

dashed line: anomalous diffusion ^' \ 

dotted line: normal diffusion /' \ 
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Figure 6: Angular profile of gamma ray intensity for the whole GCR for 10'^ SNR in 10"'' years and for 
the actual molecular gas distribution of Figure 1 ( allowed for in an approximate way ). The results of 
calculations are averaged over the galactic latitudes |6| < 0.3° and smoothed taking into account the finite 
angular resolution of the telescope of 0.1°. Normal and anomalous diffusion arc presented by dotted and 
dashed lines respectively. The results of HESS observations are also shown by the full line for comparison. 

the diffusion coefficient, in tfie absence of clear information we adopt the 'local' value 
(pertaining to the Galaxy as a whole) . 

5.2 SgrA* alone 

It is easiest to consider this region alone to start with. With only one SN, young enough 
that the particles are sufficiently confined to the region as to give gamma rays which are 
well within the point spread function ( PSF ) of the detector ~ 0.1°, the predicted flux 
can be taken from our earher work (Erlykin, Wolfendale, 2003) 

F,i> 1 TeV) ^ 35 • lO-^ (^) (^) [£j " A • / c^-.- (1) 

Here, for a start, we adopt for the SN explosion energy Eo = 10^^ erg, the gas density 
n = 100 cm~^, the distance to the GC d = 8.5 kpc and the fraction of the SN explosion 
energy transferred to GR, A = 0.1, as usual (Berezhko et al., 1996). / is the efficiency- 
factor, which is actually a ratio of the observed gamma-ray intensity to that expected 
for a standard model of the GR production, to be determined from the observations. 
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/ can, of course, be much less than unity if the physical conditions in the SNR differ 
considerably from the conventional ones, resulting in weaker shocks, etc. The result is 
F^(> ITeV) = 5/ X 10-^2 cm-^ s'K 

Following the HESS work, we adopt 0.38 TeV as the threshold energy and, using the 
measured spectrum of gamma rays (with differential exponent 7 = 2.3), the observed 
total flux from the GC including the point source at SgrA* is 
F^(> 0.38TeV) = 24 x IQ-^^ ^^-2 ^-i 
F^{> ITeV) = 7 X 10-12 ^^-2 ^-i 

Taken at face value, this would require, for one SN alone, / ~ 1.4. Thus, for 100 or 1000 
SN we would need / fa lO^^ or 10"^ It is evident that there is no shortage of energy in 
the SNR hypothesis. 

Later observations have given slightly different values of the fluxes and threshold 
energy but our arguments are unchanged. 

Figures 5 and 6 allow us to make a more accurate estimate. The flux from the spike 
within |/| < 0.2° would qualify as that from the discrete source insofar as its finite width 
would not be detected in the presence of a PSF with the half width at half maximum of 
0.1° ( which has a wide tail ). 

Of the two situations considered, it is evident that 10^ SN in 10^ years gives too steep 
a 'lateral distribution' ( comparing Figure 5 with Figure 2 ) but the shape for 10^ SN in 
10^ years is close. Here, for |/| < 0.2°, we predict a flux of 5 • lO-^cm-^s-^ for the case of 
anomalous diffusion. The observed flux is 7 • lO^^^cm^^s^^ so that / ~ 1.4 • 10"^ for this 
case. 

5.3 The Galactic Ridge 

As remarked already it is evident that the fall of predicted gamma ray intensity with 
increasing longitude is about right to explain the gamma ray profile for 10^ SN in 10^ 
years. There is a problem, however, in that the 'inner wings' of the predicted distribution, 
from \l\ : 0.05° to 0.4°, are too strong. A likely explanation here is in terms of a non- 
uniform temporal distribution. An estimate of the expected flux in the Ridge can be 
derived from Figure 6: it is 8.6 • lO-^cm-^s-^. With the efficiency factor of 1.4 • 10~^ 
derived for the sources in SgrA* we have 12 • 10~^'^ cm~^s'^, to be compared with the 
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observed 17 • 10~^^cm~^s~^. The result is close, although we have to admit that many 
parameters involved in the calculations are uncertain. 

This is where the SNR, which undoubtedly must be present in the Ridge material 
itself, may assume importance. Their likely number can be estimated, as follows. For the 
Galaxy as a whole, with molecular mass M(H2) ~ 10^ Mq ( Dame et al., 2001 ) and SN 
rate 10~^ y~^, we have ~ 10^^-*^ U^^- It is Ukely that this value is also appropriate to 
the Ridge region so that in 10^ y, with a mass of 4.4 ■ 10'^ Mq, wc expect 44. In fact, the 
number may be nearer 10 in view of the dependence of SNR density on column density 
of molecular hydrogen being slower than linear, more nearly to the power 0.6 ( using the 
summary of Fatemi and Wolfendale, 1996 for SNR, pulsars and N{H2) ). Assuming that 
these SN give the same contribution to the flux in the wings as the central SN and using 
the calculated flux in the Ridge for 1000 central SN of 8.6 • 10"^cm~^s~^ and observed 
flux of 1.7 • 10~^^cm~^s~^, the derived / - value is ~ 0.1. 

A measure of validity for the view of a significant, or even major, contribution from 
the 'conventional' SNR comes from the CR-radio correlation ( Figure 3 and §3.3 ). Three 
of the diffuse radio sources appear in the SNR catalogue of Green ( 2000 ) and in general 
we would expect the radio map to be a good indicator of past and present SNR. A similar 
situation arises for the radio map and identified SNR for the Large and Small Magellanic 
Glouds ( Mills and Turtle, 1984 ). A natural explanation would exist for the 'bumps' in 
Figure 2 in terms of SNR, in view of the expected profiles of GR-produced gamma rays 
shown in Figure 4; SNR of ( typical ) age a few 10^ years would have profiles similar 
to the 'bumps' in Figure 2. Specifically, the mean half-width of the bumps which can 
be resolved is ~ 0.24° and this is just what would be expected for about 10 SN having 
exploded in the last 10^ years. 

So far, therefore, there seem to be two possibilities to explain the results. Firstly, 
there were 10'^ SN in 10''"^ years in SgrA*, the GR diffusing through the Ridge causing the 
Ridge emission and the recent SN giving GR very close to SgrA* which caused the point 
source. The problem is that the needed gas density to get the observed ratio of Ridge 
flux to SgrA* 'point source' flux is very low, by SgrA* standards. 

Secondly, perhaps the SNR in the Ridge itself were responsible for the GR there. The 
number of SN in SgrA* could then be smaller, with higher density gas allowed. The 
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problem with the high intensity predicted in the range / : 0.05° to 3° or 4° would then be 
minimized - it would be due simply to a statistical fluctuation. 

5.4 The efficiency- factor, / 

There are many phenomena which could contribute to / being less then unity for the 
unusual conditions in the GCR. These are, mainly, 

1. The possibility that the majority of SN there are not of Type II - the main sources 
of SNR which accelerate CR to very high energies. 

2 

2. The high gas density causes the Sedov radius (which is proportional to n~5, Axford, 
1981) to be small; specifically it falls to only a few pc. The time taken to reach this 
radius, after which CR acceleration is reduced, is probably too short for efficient 
acceleration, despite the increased magnetic field in the GCR. 

3. The high gas density probably causes the injection efficiency to be low ( Drury et 
al., 1996 and private communication ), the point being that ionization losses will be 
considerable during injection for the sub-relativistic particles. 

4. The tube-like magnetic fields, referred to in §1, which will convey particles out of 
the Galaxy. 

5. A very likely effect relates to the Galactic Wind. This is currently very strong in 
the GCR (eg Breitschwerdt et al.,2002; Volk and Zirakashvili,2004; Sunyaev et al 
1993) - with a velocity of a few thousand km s^^. 

Concerning item 1, it seems unlikely that there is a shortage of the necessary massive 
pre-SN stars. Indeed, in parts (eg SgrB2), there are unusually massive stars being pro- 
duced 'furiously' ( I ). Presumably the high rate of SN production overall in the GCR, 
and particularly in SgrA*, gives rise to many Type II SN. 

The other factors are, therefore, considered to be the relevant ones. 

5.5 The Cosmic Ray Gradient 

A well known feature of the CR distribution in the Galaxy as a whole - as inferred from 
gamma ray data - is the smallness of the 'CR gradient', i.e. the dependence of CR intensity 
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on Galactocentric radius, R (eg Erlykin et al., 1996; Hunter, 2001.) and we believe that 
this has relevance, here. As is well known, the gradient for CR is far less than that for 
the (assumed) parent SNR. A number of possibilities have been put forward, involving, 
for example. Galactic winds ( Breitschwerdt et al.,2002; Volk and Zirakashvih,2004 ) and 
radial dependent diffusion coefficients ( Erlykin and Wolfendale, 2002 ). 

The results from the present analysis of the GCR lead to our suggesting that the 
efficiency of CR production by SNR compared with the standard model where 10% of the 
explosion energy is transformed into CR, is dependent on the nominal gas density of the 
ambient medium into which the SN expands. It will be apparent from the start that the 
low implied CR intensity at the GC despite the large number of SNR is in the spirit of a 
gas density dependent SNR efficiency. 

In Figure 7 we present the 'CR efficiency', given here as the ratio of inferred CR 
intensity, I(CR) to the surface density of SNR (as used by us previously, Erlykin and 
Wolfendale, 2005) versus the surface density of molecular hydrogen, a{H2), as well as the 
GC values. The fall-off of efficiency with a{H2) is very marked. 

5.6 The preferred model 

At this stage it is not possible to decide between the two models but the mixed model 
( ie SN in both SgrA* and the Ridge ) appears more likely. The /-values for the two 
situations are given in Figure 7. Consideration can be given to each region in turn. 

SgrA*. One hundred SN is probably the smallest that can be allowed, and able to 
give the necessary energy injection for the high wind, plasma temperature, etc. Thus, 
for one thousand SN the figure is 1.4 • 10~^ ( not much smaller because much of the CR 
energy escapes from the central region ). 

The Galactic Ridge. For 10 sources, we have / ~ 0.1. The values are plotted in 
Figure 7. The abscissae arc ilhistrative in that they are the approximate surface densities 
of gas for the regimes in question. Also shown in Figure 7 is the result for the source 
GO. 9-1-0.1, which appears to be due to a single SNR. It should be remarked that the region 
in which the SN exploded is outside the region of high gas density and thus the efficiency 
might be expected to be comparatively high. 
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Figure 7: Efficiency of SNR for accelerating CR of energy and intensity sufficient to give gamma rays of 
GeV energy ( CR gradient results ) and above 0.38 TeV as a function of the surface density of molecular 
hydrogen. The local values of efficiency and surface density of gas are taken as datum in each case. For 
the TeV gamma rays we adopt a nominal mean density of 100 cm~^. 



It is relevant to see to what extent the results reported here have correspondence in the 
GeV region. To this end we have examined the CGRO observations reported by Hunter 
et al. ( 1997 ); these results referring to energy ranges from 30/100 MeV to Ej > 1000 
MeV. Considering the results for |6| < 2°, the region for |/| < 1° has a spectrum flatter 
than the surroundings by A7 = 0.15, when attention is devoted to the contribution from 
interactions with molecular hydrogen alone. This is in the spirit of that reported here 
for HESS, where the spectrum of gamma rays in the GCR is flatter than average, by 
A7 ~ 0.31 ± 0.22. Bearing in mind the inferior angular resolution of CGRO and the 
quoted 'errors', the difference in the A7- values is understandable. 

Of interest, too, is the fact that there are small gamma-ray excesses for |6| < 2°, 
> 1000 MeV at values where there are also peaks in the 408 MHz maps of Haslam 
et al. ( 1981 ). Specifically, of 23 peaks in the radio map above a consistent height ( 50K 
compared with a central region of height ~450K and an Anti-Centre level of ~50K ) some 
14 coincide with gamma ray excesses. Thus, the GCR correlation of gamma-ray excesses 



5.7 CGRO data in the GeV region 
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with radio intensity seems to have a counterpart at GeV energies. 

6 Discussion 

The HESS group appear not to have analysed the results from as detailed a standpoint 
as we ourselves; instead, they take an empirical approach by way of fitting a Gaussian to 
the cosmic ray distribution ( see Figure 1 ) in GCR and a point source in the GC itself. 
This can perhaps be regarded as a zero order approach. 

Some general remarks about the analysis are in order. In §1 we discussed the general 
energetics of the GCR and pointed out the very high - and nearly equal - energy densities 
of the major components. The cosmic ray energy density is clearly lower by 3 orders of 
magnitude. Of the discrepancies, that between the magnetic field and the CR energy 
density is, at first sight, hardest to understand. Locally ( ie in the solar system ) there 
is equality and good reasons have been put forward in terms of trapping, but this must 
break down in the GCR. Presumably, the reason is the field geometry, ie the outward- 
directed nature of the high field fiux-tubes, which helps to lose CR, unlike locally, where 
the magnetic field acts to contain them. 

Another topic needing discussion is that relating to symmetry - and lack of it. In- 
spection of Figure 1 shows the distinct lack of symmetry of the density of molecular gas 
about an axis through the Galactic Centre along / = 0°. Specifically, the symmetry in the 
column density versus longitude plane ( Figure 1 ) is about / = +0.5°. This is in sharp 
contrast to that for the CR intensity, ( Figure 2 ), which is about / = 0.0° and, indeed, 
the gamma ray intensity itself ( Figure 1 ), which is about I — +0.1°. The radio intensity 
( Figure 2 ) is also symmetrical about / = 0.0°. Of particular interest is the map for 60// 
( Uchida et al.,1996 ), this radiation arising from hot dust; for this map there is symmetry 
about I — 0.0°. Most importantly, the 60// has almost disappeared by / = 1.0°, clearly, 
the dust in the anomalous region, / : 1° to 2°, is 'cold'. This result is in accord with our 
contention that the CR intensity is very low there, by virtue of the ( chance ) paucity of 
recent SN. 
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7 Conclusions 



The main results can be summarised thus. 

1. Supernova remnants can be invoked to explain the HESS results for the Galactic Cen- 
tre Region. The rate for the SgrA* region, 10^ SN in 10^ years, as recommended by other 
work ( I ), gives a consistent picture, although somewhat fewer, taken with 'conventional' 
SNR in the Ridge, is preferred. 

2. The 'consistent picture' requires particularly low SNR efficiency- factors for SgrA* ( but 
less so for the Ridge if 'conventional' SNR play a role ). Taken together with results from 
the well-known small 'CR- gradient' in the Galaxy as a whole, they lead to a satisfying 
smooth dependence of efficiency factor on local molecular density. The ideas put forward 
in §5.4 would be expected to give such a dependence. 

3. An interesting feature is the correlation of CR intensity with smoothed radio intensity. 
We regard this as providing support for the idea that SN in the Ridge are important. 

4. The CGRO data in the sub-GeV and GeV region are consistent insofar as they, too, 
give a somewhat flatter gamma ray spectral exponent for the GCR than for the rest and 
there are correlations of intensity excesses with spikes in the radio map. 

5. The very high magnetic field inhibits high energy electrons, thereby enabling conclu- 
sions about protons ( more accurately, nuclei ) to be drawn. 

Acknowledgments 

One of authors (ADE) thanks the Royal Society and the University of Durham for 
financial support. We also thank L.O'C.Drury, T.Bell, P.M.Chadwick and J.L.Osborne 
for helpful comments. Remarks by the referees of a previous version of this paper were 
particularly useful; we thank them. 

References 

Aharonian, F.A. et al., Astron. Astrophys., 425, L13 (2004) 

Aharonian, F.A. et al.. Nature, 439, 695 (2006) 

Albert J. et al., Astrophys. J., 638, LlOl (2006) 

Altenhoff, W.J. et al., Astron. Astrophys. Suppl. Ser., 35, 23 (1979) 

Axford, W.I., Proc. 17th Int. Cosmic Ray Conf. (Paris), 12, 155 (1981) 

Bania, T.M., Astrophys. J., 216, 381 (1977) 

Berezhko E.G. et al. JETP, 88, 1 (1996) 



19 



Biermann P.L., Proc. 23d Int. Cosmic Ray Conf., Calgary, Inv. Rapp. and High. Papers, 
45 (1993) 

Breitschwerdt D. et al., Astron. Astrophys., 385, 216 (2002) 

Dame, T.M., Hartmann Dap and Thaddeus, P., Astrophys. J., 547, 792 (2001) 

Drury, L.O'C. et al., Astron. Astrophys., 309, 1002 (1996) 

Erlykin, A.D. and Wolfendale, A.W., J. Phys. G: Nucl. Part. Phys., 27, 941 (2001); 
J. Phys. G: Nucl. Part. Phys. 28, 2329 (2002); J.Phys.G: Nucl. Part. Phys., 29, 641 
(2003); Proc. 29th ICRC, Pune, 3, 133 (2005) 

Erlykin, A.D. et al., Astron. Astrophys. Suppl. Ser., 120, 397 (1996) 
Erlykin, A.D., Lagutin, A.A. and Wolfendale, A.W., Astropart. Phys., 19, 351 (2003) 
Fatemi, S.J. and Wolfendale, A.W., J. Phys. G: Nucl. Part. Phys., 22, 1089 (1986) 
Green D.A., 'A catalogue of galactic supernova remnants ( 2000 August version )', Uni- 
versity of Cambridge (2000) 
Handa,T. et al., PASJ, 39, 709 (1987) 

Hartman R.C. et al., Astrophys. J. Suppl. Ser., 123, 79 (1999) 

Haslam, C.G.T., Salter, C.J. and Stoffel H., 'Origin of Cosmic Rays', ed. by G.Setti, 
G.Spada and A.W.Wolfendale, lAU, 217 (1981) 
Hunter, S.D. et al., Astrophys. J., 481, 205 (1997) 

Hunter, S.D., 'High Energy Gamma Ray Astronomy' ed. by F.A.Aharonian and H. J.Volk, 
American Institute of Physics, 1-56396-990-4, p. 171 (2001) 

Issa, M.R. and Wolfendale A.W., J. Phys. G: Nucl. Part. Phys., 7, L187 (1981) 
Kosack K. et al, Astrophys. J., 608, L97 (2004) 
Krabbe, A. et al., Astrophys. J. Lett., 447, L95 (1995) 

Lagutin A.A. et al., Nucl. Phys. B ( Proc. Suppl ), 97, 267 (2001); Proc. 27th ICRC, 
Hamburg, 5, 1900 (2001) 

Mayer-Hasselwandler H.A. et al., Astron. Astrophys., 103, 164 (1982); Astron. Astro- 
phys. 335, 161 (1998) 

Mills, B.Y. and Turtle, A.J., 'Structure of the Magellanic Clouds', eds. S.Vander Berrgh 
and K.S. de Beer, lAU Suymp., 108 (1984) 

Morris, M. and Serabyn, E., Ann. Rev. Astron. Astrophys., 34, 645 (1996) 



Oka, T. et al, astro-ph/ 98 10434 (1998) 



20 



Strong A.W. and Mattox J.R., Astron. Astrophys., 308, L21 (1996) 

Sunyaev, R., Markevitch, M. and Pavlinsky, M., Astrophys. J., 407, 606 (1993) 

Sveshnikova L.G. Astron. Astrophys., 409, 799 (2003) 

Thompson, D.J. et al., 'The Structure and Content of the Galaxy and Galactic Gamma 
Rays', p.l, GSFC, X-662-76-154 (1976) 

Tsuboi, M., Handa, T. and Ukita, N., Astrophys. J., 120, 1 (1999) 

Tsuchiya K. et al., Astrophys. J., 606, L115 (2004) 

Uchida K.I. et al., Astrophys. J., 462, 768 (1996) 

Volk H.J. and Zirakashvih V.N., Astron. Astrophys., 417, 807 (2004) 

Wolfendale, A.W., Quart. Journ. Roy. Astr. Soc, 24, 122 (1983) 



21 



